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We calculate numerically the spectrum of disordered electrons in the lowest Landau 
level at filling factor 1/5 using the self-consistent Hartree-Fock approximation for 
systems containing up to 400 flux quanta. Special attention is paid to the correct 
treatment of the q = component of the Coulomb interaction. For sufficiently 
strong disorder, the system is an insulator at this filling factor. We observe nu- 
merically a Coulomb gap in the single-particle density of states (DOS). The DOS 
agrees quantitatively with the predictions for classical point charges. 

The Coulomb interaction between charged locahzed particles leads to a van- 
ishing single particle density of states (DOS) at the Fermi energy, p{E). This 
classical effect was understood by Efros and Shklovskii as a necessary condition 
for the stability of the ground state with respect to particle-hole excitations u. A 
more quantitative description of this Coulomb gap was obtained by Efros using a 
self-consistent method □. 

The argument leading to the Coulomb gap is expected to hold in a quantum 
mechanical system, too, if the localization length of the particles is small compared 
to their distance. An important example of such a system is a disordered two- 
dimensionaLelectron gas in a strong magnetic field, exhibiting the integer quantum 
Hall effect S. In this system, all states are exponentially localized except those 
at single energies close to the centers of each Landau level Indeed, Yang and 
MacDonald numerically found a linearly vanishing DOS when the filling factor v of 
the lowest Landau level was 1/5 0. Remarkably,. the vanishing of the DOS persists 
even at the critical energy in the band center M. 

In this communication, we show that the DOS in the lowest Landau level at 
filling -factor 1/5 can be quantitatively described by the self-consistent method of 
Efros □. We carefully consider the finite-size effects in the self-consistent Hartree- 
Fock (HF) approximation and apply Efros' ideas to the finite systems considered 
in the numerical simulations. The agreement is surprising, since the effects of 
wavefunction overlap and exchange interaction are present in the distribution of 
the HF matrix elements. 

Before restricting our discussion to the QH system, we want to discuss the self- 
consistent Hartree-Fock approximation of spinless (or fully spinpolarized, as in the 
lowest Landau level) electrons on a finite torus of area L x L. The Hamiltonian is 

JV 1 ''^ 

n = J2 + t^dis(f.) + v;i_bg(f.)] + 2 E ^(^^ - + ^bg-bg, (1) 

where Ti is the kinetic energy, Udis the disorder potential, V^i-bg the interaction 
with the homogeneous background charge, V the Coulomb interaction, and Vbg^bg 
the background self-interaction, respectively. The Hartree-Fock equations then take 
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the form 

(r + {/dis + Ki-bg)^a(r) + 5]//3 5] [(/3|/3)q^a(r) - (/3|a)qV/3(r)]z;,e'9r = sM^), 

P q 

(2) 

with (/3|a)q = {f3\e *'i ''|a), Vq — 2TTe'^ / {L^q) is the two-dimensional Fourier trans- 
form of the Coulomb interaction, and fa is the Fermi function. The q = part of 
the sum over q gives 

V0^fl3'fpa{r){l - 6a, 0) Vo{N ~ fa)i}a{r) = - V^l-bg^Aa (r) - Wo/aVa(r), (3) 

which cancels the interaction with the background charge up to a contribution 
proportional to vq — JcPrr^^L^^ ^ L^^ < oo. Since this contribution vanishes in 
the limit of infinite system size it is often neglected resulting in the HF equations 

(T + [/dis)^a(r)+ J2 fp[{mMr)-ma)^i^p{r)]vge''^^ ^EaMr)- (4) 
However, the correct eigenvalues for the finite system would be 

Si^ Ei- fiVQ, (5) 

as an electron in an empty state interacts with N other electrons, while an electron 
in an occupied state interacts with only N j^l other electrons. 

Efros based his self-consistent method □ on the observation that the ground 
state of a system is stable against every particle-hole excitation. For the solutions 
£i of eq. this is equivalent to the condition 

AE = Ej - £, - > 0, (6) 

for every In the classical limit, for states localized in an infinitesimal vicinity 
of the position r^, the interaction matrix element is given by 

«»i ^ZZ^I («l*)q(ib')-q = I 1^ I - ^^0- (7) 

For the torus geometry, 

vo = Y^ dV- =21n(l + ^)- « 1.76^. (8) 
^ JlxL t L L 

In the classical limit, eq. restricts the possible distances jr^ — Yj \ and energies 
£j for a given energy e^. This leads to a self-consistency condition for the DOS □ 



p(e) = po exp 



pje') de' 
2" Jo {e + e' + voy 



(9) 



In contrast to the infinite system result of Efros, eq. in a finite system contains 
the finite contribution vq- 

In the limit of small e and large L, eq. has the asymptotic solution 

p{e) = -^{vo + \s\). (10) 
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The stability condition for the true HF eigenvalues ii (eq. (||)) is 

/S.E — Sj — £1 — Vij — vq> 0. (11) 

The self-consistency equation for the density of the e does not contain wq. However, 
now a lower bound for the energy difference ij — £i exist, since the distance between 
the states can be at most L/\/2 on the torus: 



p{£')d£' 
JV2eyL {£ + £'? 



(12) 



p{e) = (aV2^ + |£| ) , (13) 



p(e) = pq exp 
Again, the asymptotic solution is 

where a « 0.567 is the solution of In a = —a. 

We solved the self- consistent HF equations for lowest Landau level with a short- 
ranged disorder potential u u. The resulting DOS for a system containing 400 flux 
quanta at filling factor 1/5 averaged over 60 disorder realizations is shown in Fig. 0. 
In order to avoid a filling in of the gap due to fluctuations of the Fermi energy, the 
spectra were shifted to assure Eat -f-ejv+i — for every disorder reahzation, where N 
is the number of electrons. Comparison with the self-consistent solution of eq. (^ 
and its asymptotic solution eq. (|lO|) shows qualitative agreement. The deviations 
at larger distance from the Fermi energy are due to the assumption of a constant 
DOS in the derivation of eq. (^. 

The DOS of the true HF eigenvalues e shows qualitatively the same behavior. 
In Fig. |2| a) the DOS at the Fermi energy is shown for different system sizes and 
disorder strength. The data depend only on lB/{"fL), where 7 — e'^/{lB^) is a 
measure of the strength of the interactions relative to the disorder F, and for small 
arguments are well described by the asymptotic solution eq. (|l3|). 

The agreement between numerical simulations and the classical self-consistent 
method is surprising, as the conditions for its applicability are only rather poorly 
fulfilled in our system. The localization length ^ m 41b oi the non- interacting system 
is only slightly smaller than the mean particle spacing of (47r/(3^/^i/))"'^/^ZB ~ 61b 
Thus one might expect to see stronger effects of the finite wavefunction overlap and 
the exchange interaction. Indeed, the distribution of interaction matrix elements 
shows these effects. Corresponding to a matrix element Vij we define an effective 
distance r^B = ^ Ik^ij + ^0) and plot its histogram in Fig. |^ b). For a random 
distribution of classical point charges one would expect a distribution that rises 
linearly for small distances and that terminates at the maximum distance L/-\/2. 
In contrast, effective distances much larger than L are present corresponding to 
states with small interaction matrix elements due to the exchange interaction. 

We studied numerically the Coulomb gap at low filling factor in an interacting, 
disordered model of the integer quantum Hall effect. Surprisingly, we find quan- 
titative agreement with the predictions of Efros' classical self-consistent method, 
despite the fact that the distribution of interaction matrix elements shows strong 
effects of the quantum mechanical effects of wavefunction overlap and exchange. 

This work was performed within the research program of the SFB 341 of the 
DFG. 
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Figure 1. The density of states obtained by smearing every eigenvalue by a Gaussian of width 
0.8 mean level spacings for interaction strengths e^/(isr) of 0.025, 0.15, 0.275, 0.4, 11.525, and 
0.65 (top left to right bottom). The dashed lines correspond to the solution of eq. W) and the 
asymptotic solution eq. dlQ). 




Figure 2. (a) The averaged DOS at the Fermi energy p(0) = l/(Z/^(ejv+i — ejv)) of the true HF 
eigenvalues as a function of lg/{'yL). (b) Distribution of effective distance r^g of the HF states. 
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